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Abstract It is shown that the description of collapse given by the classic 
model of Oppenheimer and Snyder fails to satisfy a crucial matching condition 
at the surface of the ball. After correcting the model so that the interior and 
, exterior metrics match correctly, it is established that the contraction process 

stops at the Schwarzschild radius, that there is an accumulation of particles at 
the surface of the ball, and that in the limit of infinite time lapse the density of 
particles at the surface becomes infinite. A black hole cannot form. This result 
confirms the judgements of both Einstein and Eddington about gravitational 
collapse when the collapse velocity approaches that of light. 

1 Introduction 

The Relativistic Theory of Gravitation (RTG)[T] may be considered an exam- 
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pie of a class of theories introduced by Rosenp], known as bimetric theories. 
Coexisting with the field, or Ricmann metric 

ds 2 = g^dx^dx" , (1) 



J> ' there is a space, or Minkowski metric 

OS " 

CO; dn*=^dx*dx 1 ' . (2) 

Unlike in General Relativity (GR), where the field is synonymous with the ge- 
C ometry, we recognize that the field, as in the electromagnetic case, is propagated 

through an underlying Minkowski space. In RTG there is a preferred coordinate 
system, namely that for which the space metric is Galilean. It may be called 
the inertial system or frame, and its corresponding field metric g^ v gives rise to 
a gravitational potential = g^yf—g whose divergence is zero. The latter 
condition may be made covariant by requiring the covariant divergence of the 
gravitational field in the Minkowski metric to be zero. This is the coordinate 
system which Einstein [3j used in order to derive his formula for the gravitational 
radiation emitted from a time-varying quadrupole source. 

In this article I shall show how the comoving coordinate frame introduced 
by Tolman[4], and developed by Oppenheimer and Snyder[5][6] (OS) to describe 
gravitational collapse, may be transformed to the inertial, or harmonic frame. 
It is then possible to track the trajectories of individual particles in the collapse 
of a dust ball, for which the equation of state is simply p = 0. We establish 
that the collapse can go no further than the Schwarzschild radius, that the col- 
lapse to this state takes an infinite time, and that the density of particles in the 
limit becomes infinite at the surface of the ball. This confirms the description 
given by Logunovp. , in which gravity changes from being attractive to repulsive 
for certain high-density conditions. The description of OS, namely "continued 
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gravitational contraction", which may be considered the ancestor of the con- 
temporary "black hole", is now seen, in the light of RTG, to be incorrect. It is 
more appropriate to call the process one of gravitational compression; the com- 
bination of an attraction of the surface particles with a repulsion of the particles 
beneath the surface produces an infinite density at the surface; presumably this 
balance of purely gravitational forces will be modified by the action of nuclear 
forces in an actual neutron star, so the state of infinite density is never reached. 
But, in any case, contraction does not go beyond the Schwarzschild radius, as I 
shall be able to show by tracking the trajectories of individual particles. The re- 
sults obtained here confirm the intuitions of both Eddington 7 and EinsteinJS], 
as expressed in the early period of the black-hole hypothesis. Particles which 
crossed the "event horizon" would acquire velocities exceeding that of light, 
which runs contrary to the prescription of Einstein's Special Theory of 1905; so 
black holes are ruled out by that theory. 

It may be the case that a modified Einstcin-Hilbcrt equation incorporating 
a nonzero cosmological constant, a natural and possibly even necessary fea- 
ture of RTG, will, as suggested in a recent article of Gerstein, Logunov and 
Mestvirishvili[9], stop the compression process after a certain time and convert 
it into a cyclical, or pulsating, process. However, it is not clear what will happen 
if, as seems plausible, nuclear forces intervene before these extra gravitational 
terms can take effect. 

Although RTG is a distinct theory, constructed in opposition to GR and 
following ideas planted by earlier workers like Rosen [2] and Fock[10], it is pos- 
sible to view the results obtained in the present article as simply a way of 
transforming a long established solution of the GR field equations, namely that 
of Tolman, into a more realistic coordinate system. Previous attempts of this 
nature, starting from the Schwarzschild solution, are those of Eddington and 
Finkelstein and of Kruskal, described in the book of Landau and Lifshitz[6] (see 
also HartlepJ]), but these made no attempt to penetrate into the interior of the 
collapsing object. 



2 The Oppenheimer-Snyder metric 

The Oppenheimer-Snyder (OS) field metric for a dust ball may be written 

ds 2 = air 2 - V 2 dR 2 - W 2 (dO 2 + sin 2 0d(f> 2 ) , (3) 



where 

2/3 



2m \ 4m '/ 2m ^2mW 



(4) 



In this comoving metric the free-fall radial geodesies are simply R =constant, 
and the coordinate r is the particle's proper time; in particular the surface of 
the ball is specified by a fixed R, which will be put as 1. The function F (R), or 
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rather the product FF' , gives the mass distribution of the dust particles. For 
R > 1 (the external region), we put 



F(R) = 1 (R>1) 



(5) 



giving zero density there, and for R < 1 (the internal region), F is left arbitrary 
for the moment, but we note that, for the metric to be continuous it has to 
satisfy 

F(l-) = 1, F'(l-) = . (6) 

At this point OS made a fatal error by choosing an F which fails to satisfy the 
second of these, so my choice of F will constitute a corrected version of OS. The 
cumulative mass distribution is given by the function 



M{R) 



AttT 00 (R') y/-g(R')dR' = mF 2 (R) 



(7) 



The harmonic coordinates [10] [12] [M] (t,r,9,<p) for this system are given by 
the solutions of 

□t = 0, Dr=~ , (8) 
where the spherical d'Alembertian operator is given by 
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and we use dot and prime to signify differentiation with respect to r and R 
respectively. The solution will be chosen to satisfy the asymptotic conditions 



t ~ r, r ~ W (r 



-oo) 



(10) 



and in that case the evolution R =constant describes, in the asymptotic re- 
gion, the Newtonian collapse of a dust ball of uniform density. This may be 
demonstrated from the relation between r and t, which is 



9m\t\ 2 F 2 (R) 



r(R,t) 

which, combined with ([7]), gives 

M(R) = m 



1/3 



(i -oo) 



r(R,t) 



•(1,*) 



(11) 



(12) 



and shows that the mass contained in a ball of radius r is proportional to the 
ball's volume. The t-dependence of r has, of course, been known since Michell 
discovered it in the eighteenth century. 

In the exterior region the solutions of JH]) are[14j 



t = T- 2\j2mW + 2m In 



W 
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'2m 



r = W-m (R>0) 



(13) 
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and the OS metric becomes, in this region, 

dg2 = r-ffl^j _ r_+jn dr2 _ 2 ^2 + sin 2 ^ (M) 

r + m r — m 

which is the Schwarzschild metric, except that the Schwarzschild "radius" r has 
been replaced by r + m. We shall see that the solution of ((8]) may be continued 
into the interior region all the way to R = 0, and that this corresponds to r = 
for all t. I deduce that it is this r, rather than Schwarzschild's, which should 
be regarded as the true radius, a conclusion which will be reinforced by making 
a closer examination of the internal solutions of ([8|). It is clear that the above 
solution (fT3|) gives, as t goes to plus infinity, that W approaches 2m and r 
approaches m, and that these limits are reached at a finite value of proper time 

Tf (R) = ^(VB?-l) (R>1) . (15) 

The fact that this is finite has led to the widespread, indeed almost universally 
held, conclusion that a falling particle goes on to cross the "event horizon" at 
r — m for proper times greater than t/, and therefore to be swallowed by the 
black hole at r = 0. I shall show that this conclusion is incorrect. 

All I have to do is demonstrate that the interior solutions of © give a 
limiting value r/ (R) in < R < 1, with rj (0) = and 77 (1) = m, together 
with a corresponding t f (R) , for which t goes to plus infinity. I have been able 
to make the demonstration numerically for rj (R) with a particular choice of the 
function F (R) . However, for general F, the values of r/ (R) may be obtained 
simply by examining the characteristics of © . These satisfy, for both of these 
partial differential equations, the same pair of ordinary differential equations 

dr / / — ,\ /2m 

— =±V(T,R)=±(2mVR-TF') S J— . (16) 

It is a simple matter to verify that the above value of t f (R) satisfies this with 
the upper sign in R > 1, where F = 1 and F' = 0. In the interior region the 
upper-sign characteristic through the point 

=0 (17) 

is the first one to be met by a geodesic R =constant coming from r = — 00. The 
situation is illustrated in Figure 1, where I have plotted this characteristic, in 
both exterior and interior regions, for the case 2m — 1, using in the interior the 
function F of my previous article [T5], namely 

F (R) = R 3/2 e 3X/2 , X=l-R (R < 1) . (18) 

The figure represents the whole of space, that is < r < +00, and the whole 
of time, that is -00 < t < +00; there are no singularities and no trapped surfaces 
of the kind supposed by Penrose and Hawking |15j. Such strange objects are to 
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Figure 1: The limit of physical space-time with 2m = 1. r is the proper time 
of a dust particle and R is its comoving coordinate, so that R = 1 indicates a 
surface particle. A given dust particle, or in the exterior region a test particle, 
moves along the abcissa R =constant, arriving at the boundary curve after an 
infinite time t. 
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be found in the region beyond t = +00, and belong in the realm of science 
fiction. The finite value of the proper time T/ is an indication that a falling 
particle, as it approaches the boundary characteristic passing through the point 
R = 1, r = 0, suffers an infinite gravitational red shift. We shall see a natural 
explanation for this in the infinite surface density of the dust ball in the limit 
t — > +00. 

The integration of ([8]) in the interior region may be best achieved by changing 
to the characteristic coordinates (R, S) , defined by 

t = t(R,S), T(l,S) = -2mS, Drt = V = UmVR-TF^j 

(19) 

The family of curves given by S =constant are the upper characteristics of 
©, and in particular the boundary characteristic is S = 0. In terms of these 
coordinates the d'Alembertian ^ is transformed by putting 

D R = d R + Vd T , (20) 

where D R is the derivative with respect to R along the characteristic, giving 

1 o 2 f fV 2V\ 2 / • \ 

V = -y,Dl+ V drD R + v ,( y - w )D R+w (l + w)d T . (21) 

Note that d T does not commute with D R , but that if we write it as 

d T = *- l D s , *=d S T(R,S) , (22) 

then D R and Ds commute. In the numerical procedure, described in the next 
Section, there is no need to find the function \P explicitly, because we work with 
the differential operator d T rather than Ds- Although this integration is rather 
formidable, it may be guaranteed that no singularity occurs in the whole physical 
region, because the coefficients of the PDEs are nonzero and nonsingular, and 
the determinant of the second-order coefficients is negative throughout, thereby 
preserving their hyperbolic character. 

Note that the characteristic curves we just introduced are also the null 
geodesies of the (corrected) OS metric. They are the paths taken by outward 
going light signals, and since the local light velocity is V, their form underlines 
the importance of the correction I made to the original OS metric, insisting on 
V being continuous at R = 1. The time required for such a signal to go from 
an internal point (Ri, S) to an external point (R2, S) may be written as 

t(R 2 ,S)-t(R 1 ,S) = t(l,S)-t(R 1 ,S) + G(R 2 ,S) , (23) 

where G is the travelling time in the exterior region, that is 

G(R 2 ,S) = t(R 2 ,S)-t(l,S) . (24) 

In view of the simple form (|14[) of the exterior metric, this latter integral may 
be simplified to give 

G (R 2 , S) = r (R 2 , S)-r (1, S) + 2m In - ^f 2 ' ^ - m . (25) 
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Now, since r(l, S) > m for all t, this result shows that any light signal emitted 
from inside the ball eventually reaches the exterior region, and this should lay 
to rest all preexisting ideas regarding trapped surfaces. But note that I said 
"eventually"; owing to the infinite red shift suffered at the boundary R = 1, 
this travel time becomes infinite as r (1, S) approaches m, that is in the closing 
stages of the compression process. 

3 Numerical integration 

It is convenient to define the operator 

P = -V 2 D = D 2 R - 2£Wd T D R - 2?W (l + W) d T + ^ - |) D R , (26) 
where £ = V/W, that is, for the choice we made in (|18p for F, 

£ = R- 1 (l + -Xe 3X / 2 ) ( 1 - ^ e 3X / 2 ) 1 , (27) 
\ Am J \ 4m / 

and the harmonic coordinates then satisfy 

Pt = 0, (P-2i 2 )r = . (28) 

These PDEs must be integrated inO<i?<l,S f >0 with surface boundary 
conditions at R = 1 

t (1, S) = -2mS - 2 v / 2mW Q + In ^1 + r(l,S) = W -m , (29) 

VWo — V2m 

where 

/ oc\2/3 

Wo(S) = W{l,S) = 2m\l + — \ , (30) 

and also with the asymptotic condition (|10[) for S — > oo and with finite values 
at i? = 0. 

Because t(l,S') becomes infinite at S = 0, the integration must be over a 
range 5 > Si > 0, and we may then convert each PDE into a set of coupled 
ODEs in R, starting from t(l,S) and r (1, S), for a discrete set of N values of 
S between Si and an upper limit S2 for which the asymptotic values may be 
used. Because S is itself defined by the first-order ODE p^|) . we retain r as 
the dependent variable, so there are a total of 3./V coupled ODEs for t, Dt and 
r. Fuller details, including the asymptotic matching procedure, are given in the 
Appendix. 

From the solutions r (R, S) and t(R,S) I interpolated to obtain r(R,t). 
Putting 7'o (t) — r (1, t), the relative position of a dust particle whose position 
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Figure 2: Evolution of the particle distribution u(R). u is the position of a 
given particle relative to the surface, indexed by its comoving coordinate R. 
The lower bold curve gives u{R) in the limit r — > oo, that is t — > — oo, while 
the upper bold curve gives u (i?) for ro = m, that is t — > +00. The lighter curves 
give u(R) for the values (i) r = 2m (ii) r = 1.4m (in) r = 1.1m. 

in the ball is indexed by (R, 9, <p), with < R < 1, may then be described by 
the coordinate 

u(R,r ) = r[R ' t{r ° )] (0<R<1), . (31) 

I have plotted, in Figure 2, u(R,ro) against R for various values of ro. In the 
early stage of collapse the slopes at both R = and 1 increase as ro decreases, 
indicating that particles near the centre move towards the surface and particles 
near the surface move towards the centre. However, towards the end of the 
process the curve u(R,r ) almost immediately crosses the curve u(R, 00) (the 
lower bold curve in Figure 2) near R = 1, and in the limit r — > m approaches 
the upper bold curve. The latter has both zero slope and zero curvature at 
R = 1, indicating that the particle distribution in the final state has an infinite 
density at the surface. 
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4 Conclusion 



A surprising feature of Figure 2 is that the density of dust particles, which 
started off uniform in the Newtonian region of t, becomes infinite at the sur- 
face of the ball as t — > +00; this may be considered as confirming Eddington's 
intuition [7] that something intervenes to prevent the "absurdity" (to quote Ed- 
dington) of the ball collapsing to a point. We see that the finite radius of the 
end state is a consequence of a purely gravitational field, without the need for 
any other forces to counter the gravitational attraction. Indeed, if we commit 
the heresy of speaking of gravitational forces, we may say that they turn from 
being attractive to repulsive in the high-density region; it is the accumulation of 
gravitational energy inside the ball which prevents it from collapsing to a point. 

Such heresy is at least hinted at by the favoured status afforded here to the 
harmonic coordinate frame, often referred to as a "gauge". But is it possible 
to put the results I have established into a coordinate-free form? My central 
result is the form of the boundary characteristic given by (fTI)]) and (jTTJ). This 
makes no reference to the harmonic frame, even though I went via that frame 
to establish it, so a coordinate-free derivation may be possible. 

I would, nevertheless, argue that use of the harmonic frame advances our 
understanding of the collapse, and especially of the stable final state, which 
was declared nonexistent in the OS analysis. There was a mystery about the 
infinite red shift suffered by a free-falling exterior particle as it approaches r = 
m. Now, in the light of the discovery that particles inside the ball behave in 
a similar manner, such behaviour becomes less mysterious. We merely have 
to find an explanation for the build up of density at the surface, and that, 
expressed in suitable coordinates, is to be found in the Einstein-Hilbert field 
equations themselves. It may be that some change of frame could transform 
away such an infinite density so that it becomes finite, but the theory of gravity 
which motivated the analysis of this article [1], and also the earlier statement of 
opposition to Strong Equivalence by Fock[10], seem to offer substantial support 
to the assertion that the harmonic, or inertial, coordinate frame should have a 
privileged status. 

Finally I remark that the harmonic frame was first accorded such a status in 
Einstein's pioneer article on gravitational waves[3]. The subsequent history [Tfi] 
of gravitational waves is a process of repeated assertion and denial of the Strong 
Equivalence Principle that all frames are equivalent, beginning perhaps with Ed- 
dington's observation [T7] that, if we allow arbitrary coordinate transformations, 
the waves may be made to appear and disappear at will; gravitational waves 
travel at the speed of thought. The converse of this statement, of course, is that 
if we do not allow such transformations, but instead insist on a privileged frame, 
such waves do exist and are worth searching for. 
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5 Appendix 



The values of t and r in the asymptotic region may be obtained in terms of the 
variables (R, v), where 

u 3 = e -3x/ 2 _3r x = l _ R (32) 
Then the operator P, defined by (p?fi|) . becomes 



* 2 +' ^j.-^mL -J , (33 ) 



i?X 2w 3 / i?v 3 £ \xJ V 2i? 
where 

£=§ + Jp V = er^ • (34) 

By expanding this in inverse powers of v and substituting in |28|) , one obtains 
the first few terms in the asymptotic series of t, with the first term given by the 
Newtonian limit, that is 

_^^ + J2t n (R)v 1 - n , (35) 



n=0 



with 



to = --R 2 e 2X --, tl = -e- 3X / 2 , 
2 2' 3 

. _ 91 1 E>2 2X 1 p4 AX 

f 3 = ±e~ 3X / 2 + ^(U-6X + X 2 )e x / 2 -20 , 

571 91 p2 2X , 1 p 4 « , 1 r>6 e 6X 

u vi2-m Re + im Re + 3mo Re ' 

*6 = g* 3 + 120 ( e_3X/2 - : ) + * 51 - t52 + ~R~ ib3 + t52 (0) ~ t53 (0) ' 

te = ^ + ^(9-4X + X 2 ) e -^-fe- 3 ^ 2 + i e -^ . (36) 



where 



_ 1 Xl7 2 

<51 -6i IT 



and 



1 

5 



l 

r]U 5 . (38) 
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Note that the latter three integrals may all be given in terms of elementary 
functions, though they are rather complicated. 
The corresponding series for r is 



r 
2m 



= e x R 



n=0 



n (R) V- 



(39) 



with the coefficients 



ro 


= -R 2 e 2X 
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3 

- |, r x =r 2 =r A 


= f5 = r-j = 


1 




= 20 ~r 


-3X/2 - | (13 - 6X 


fl 2 )e x ' 2 , 






5 15 

6 ~ T ' 


9-4X + X 2 )e~ x - 


h 3 


JL-3X 

12 


rs 


2403 
s 


- ^e* (X 4 - 9X 3 4 


- 40X 2 - 89X - 


1-89) 



rg 



X 2 + X + 1 



16 

675 + 
400 



20 c -3X/2 



9 96 
6^ 2 - —X 



l) 2 + 60e x/2 (X 2 - 5X + 10) 
1158 s 



25 



-5X/2 



-3X 



10 
27 f 



-9X/2 



(40) 



In each coefficient of both expansions two constants of integration had to be 
specified. They were determined from the continuity condition at R = 1 and 
the fmiteness condition at R — 0. The derivatives of t n at R — 1 are fixed by 
this procedure, and we find thalQ 



t' —t'—t' —t'—t'— n t' 

c — E 2 — l 3 — l 4 — l 6 — u ' c l 
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(41) 



Rather remarkably, although the asymptotic series breaks down for X — O (v 3 ) , 
a match with the solution in this boundary layer confirms the above values, and 
gives the additional ones 



68 

4 = 0, t' 7 = —, t' s = -0.0616. 



4 = 1 



*i Q = 0.6314. 



t' n = -2.3253 . 



(42) 



The corresponding procedure for the r coefficients gives r' n = for n < 9, and 
the boundary layer matching not only confirms these, but shows that r[ , r' n , r' 12 
are also zero. I have been emboldened to conjecture that the partial derivative 
of r with respect to R is zero for all v > 1. 



x Note that this result is at variance with the assumption of continuity for tn made by me 
in Ref|13|. Essentially, we must replace this condition by one of fmiteness at R = 0. But I 
stress that the condition imposed by Ref[5] on the interior metric in the entire region R < 
cannot be satisfied by any solution with an acceptable degree of continuity at the surface. 
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The numerical integration procedure I have described requires a knowledge 
of not only t and r at R — 1, but also their i?-derivatives there. For r the 
conjecture I just described may be tested by moving the lower limit of the region 
of integration close to S = 0, and when this was done there was no noticeable 
tendency for the solution to diverge as the limit R = was approached; of course 
rounding errors always allow the divergent solution to enter at some stage, but 
this tendency seems to be constant over the whole range of S. In the case of 
t, all we can do is include all the asymptotic coefficients we have, in order to 
approximate the initial value of its i?-derivative. Such a procedure gives a stable 
solution down to a dustball radius ro of around 1.1m, and thanks to our rather 
precise r-solution we already know what happens for infinite t, that is = m. 
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